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ABSTRACT 

We present a morphological analysis of 17 X-ray selected clusters at z ~ 0.25, imaged uniformly with Hubble 
Space Telescope WFPC2. Eight of these clusters comprise a subsample selected for their low X-ray luminosities 
10'*'*ergs^'), called the Low-Lx sample. The remaining nine clusters comprise a High-L^ subsample with 
Lx > 10"^^ ergs s~'. The two subsamples differ in their mean X-ray luminosity by a factor of 30, and span a 
range of more than 300. The clusters cover a relatively small range in redshift (z =0.17-0.3, (j,/z ^ 0.15) and 
the data are homogeneous in terms of depth, resolution (0."17= lh^^^ kpc at z = 0.25) and rest wavelength ob- 
served, minimizing differential corrections from cluster to cluster We fit the two dimensional surface brightness 
profiles of galaxies down to very faint absolute magnitudes: M702 < — 18.2 + 51og/i5() (roughly O.OIL^) with 
parametric models, and quantify their morphologies using the fractional bulge luminosity (B/T). Within a single 
WFPC2 image, covering a field of 3' (ih^^ Mpc at z = 0.25) in the cluster centre, we find that the Low-Lx 
clusters are dominated by galaxies with low B/T 0), while the High-Lx clusters are dominated by galaxies 
with intermediate B/T 0.4). We test whether this difference could arise from a universal morphology-density 
relation due to differences in the typical galaxy densities in the two samples. We find that small differences in 
the B/T distributions of the two samples persist with marginal statistical significance (98% confidence based on 
a binned test) even when we restrict the comparison to galaxies in environments with similar projected local 
galaxy densities. A related difference (also of low statistical significance) is seen between the bulge luminosity 
functions of the two cluster samples, while no difference is seen between the disk luminosity functions. From 
the correlations between these quantities, we argue that the global environment affects the population of bulges, 
over and above trends seen with local density. On the basis of this result we conclude that the destruction of disks 
through ram pressure stripping or harassment is not solely responsible for the morphology-density relation, and 
that bulge formation is less efficient in low mass clusters, perhaps reflecting a less rich merger history. 

Subject headings: galaxies: clusters: general - galaxies: structure - galaxies: evolution 



1. INTRODUCTION 



The local and large-scale environment is known to be an im- 
portant determin ant of many galaxy properties, in partic ular star 
formation r ates ([Osterbrock 1960|; [Dressier et al. 1985t B alogh 
et al. 1998; |Moss and Whittle 2000D, gas content (Bahcall 



1977; Piovanelli and Haynes 1985|; [White and Sarazin 1991 
Bravo-Alfaro et al. 2000|; [Vollmer et al. 200U Solanes et al 



2001) and especially morphology (Rubble and Humason 1931 
Spitzer and Baade 1951[; [Vlorgan 196U [Abefl 1965[; D ressier 



1980; postman and Geller 1984 ). These studies all suggest that 
galaxy morphology, in particular the presence of a star-forming 
dis k, is a fundamental defining characteristic of galaxies (Hub- 
ble [l92|;|l92|). 

The variation in morphological mix with local galaxy den- 
sity was first quantified by the work of Dressier (1980), who 
showed that the proportion of early-type (elliptical and lenticu- 
lar) galaxies increases at the expense of late-type spiral galax- 
ies in regions with higher projected galaxy density. This trend 
is called the morphology-density relation, hereafter T-E, and 
in this paper we investigate whether the T-S can provide any 



insight into the origin and permanence of galaxy morphology. 

The advent of high-resolution imaging from the Hubble 
Space Telescope (HST) has allowed greatly improved quanti- 
tative analyses of galaxy morphologies, especially beyond the 
very local universe (z > 0. 1). This has resulted in a renewed in- 
terest in the use of morphology as an important tracer of galaxy 



evolution (e. g. Co uch et al. 1994; Abraham et al. 1996a ; 
Small et al. 1997[ ; Schade, Barrientos & Lopez-Cruz [1997| ). 
For example, working from an extens ive H ST survey of rich 
clusters at z ^0.4-0.5, Dressier et al. (1997) have shown that 
the T-S seen in local clusters is present in the most massive 
systems at z 0.5. However, these studies have also suggested 
that a galaxy's morphology may evolve, complicating the com- 
parison of populations at different epochs. 

Theoretically, there are several physical mechanisms which 
may determine or alter a galaxy's morphology as a function of 
l ocal density and/or cosm i c time. Early numerical sim ulations 
( Roos and Norman 1979|; [Farouki and Shapiro 1982[; B arne s 



1992) demonstrated that a merger between two spiral galaxies 
would almost always result in a galaxy with structural parame- 
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ters resembling those of elliptical galaxies. Thus, environments 
rich in merger events (dense regions with low velocity disper- 
sion) are expected to be fertile regions for the formation of ellip- 
tical galaxies. This is an important ingredient in mo dern models 
of galaxy formation (|Cole 1991 ; [Kauffmann 1996|; Kauffmann 
et al. 1993; paugh et al. 1996| ; [Somerville and Primack 1999|) 



which produce a rea sonably good match to the observed T-S 
(Diaferio et al. 2001 ). However, there are other physical mech- 
anisms, generally not included in these m odels, which might 
affect galaxy morphology. Moore et al. ( 1999[ ) showed that 
tidal interactions between galaxies in dense environments are 
effective at destroying galactic disks and at transforming spiral 
galaxies into the dwarf spheroidal galaxies that dominate rich 
clusters today. Other p rocesses, such as ram pressure stripping 
( Gunn and Gott 1972|) turbulent and viscous stripping (Nulsen 



1982) and strangulation ( [Larson et al. 1980| ; [Balogh et al. 206o| ) 
result in the fading of galactic disks and in the smoothing out of 
their surface brightness due to the fading of HII regions. These 
processes, then, migh t transform spiral ga laxies into SO galax- 
ies or anemic spirals (van den Bergh 1991). 

One of the stronger trends seen in the T-E at low redshift is 
the dramatic increase in the number of lenticul ar (SO) galaxie s 
in the highest density regions of the universe (Dressier 1980). 
There is now also evidence that the fraction of SO galaxies in 
rich c lusters evolve s dram atically with redshift (Dressier et al. 



1997; Fasano et al. 200C, but see also Andreon 1998 and Fab 



ricant, Franx & van Dokkum 200C ), suggesting that the domi- 
nant process in defining the T-E is the transformation of spiral 
galaxies in to SO galaxies within high density regions (P oggianti 
et al. 1999; p<-odama and Small 2001 ). This suggestion appears 
to be consistent with studies of age differences between ellip 



tical and SO galaxies in clusters (Kuntschner and Davies 1998 
van Dokkum et al. 1998| ; [Jones et al. 200q ; [Small et al. 2001[ ) 



However, there remain significant questions about the feasibil 
ity of this transformation. Most importantly, it has been sug- 
gested that the correlation between bulge luminosity and bulge- 
to-disk ratio for spiral galaxies means that it is difficult to form 
the entire cluster SO population by stripping the disks of spi- 
rals, as there is too little lumino s ity density in spiral bulges 
( Yoshizawa and Wakamatsu 1975[; [Dressier 1980[ S olanes and 
Salvador-Sole 1992; [Kodama and Small 2001[ ). Since the ob- 
served evolution in this galaxy population occurs at relatively 
modest redshifts (z < 0.5), this is an issue which can be ad- 
dressed with present observational facilities. 

Another unanswered question is which physical processes, 
acting in which environments, are responsible for establ ishing 
the morphological mix of a galaxy population. Dressier ( [198C ) 
showed that galaxy morphology is strongly correlated with pro- 
jected local galaxy density, suggesting that morphology must 
change as the local density changes, irrespective of whether or 
not the galaxy is in a rich cluster However, this has been dif- 
ficult to establish with certainty because, particularly in rich 
clusters, local galaxy density is closely correlated with cluster- 
centric distance. It has been suggested that, with an appropriate 
choice of the cluster centre, morphology correlates better with 
radius than with local density and hence this global property is 
more likely to be driving the variation ( Whitmore and Gilmore 
1991; [Whitmoreetal. 1993^ 

To address this question observationally, it is critical to inves- 
tigate how galaxy properties depend on global cluster proper- 
ties such as mass or X-ray luminosity, independent of any other 
factors, such as local density or redshift. One can then hope to 



distinguish between physical mechanisms which affect galaxy 
properties on very local scales from those which are coupled to 
the la rge-scale environmen t. For example, ram pressure strip- 
ping ( |Gunn and Gott 1972 ) is strongly dependent on the galaxy 
velocity in the presence of a dense intracluster medium, and is 
thus more likely to occur within massive clusters with high ve- 
locity dispersions. On the other hand, galaxy merging is most 
effective in regions of low velocity, so may be more important 
in groups and low mass clusters. Finally, disk galaxies which 
end up in massive clusters may have intrinsically larger bulge- 
to-disk ratios than the general field population. This could oc- 
cur since it is expected that galaxies in regions which are over- 
dens e on large scales form earlier than average (B ardeen et al. 
1986; [Narayanan et al. 200C| ). As a result, such galaxies may 
have consumed more of their gas and be more likely to have re- 
duced or ceased star formation in their disks by the present day; 
consequently, the galaxies would be more bulge-dominated. 

To address these issues, we have obtained Hubble Space 
Telescope (HST) WFPC2 F702W (Rjoi-hand) imaging of 17 
clusters at 0.17 < z < 0.30 (with a mean redshift of z = 0.25), 
selected from two published ROSAT cluster samples. The 
choice of redshift is a compromise between the need for both 
deep and wide images (to ensure coverage of a large physical 
radius); furthermore, z — 0.25 represents an epoch where a sig- 
nificant SO galaxy population is visible in the cores of rich clus- 
ters, and must have been very recently transformed. Therefore, 
observations at this epoch represent the best opportunity of un- 
covering the dominant environment- related processes re sponsi- 
ble for th e creation of SO galaxies ([Jones et al. 200C|; S mall 
etal. 2001). 

All clusters were selected on the basis of their X-ray emis- 
sion. As part of a more extensive imaging and spectroscopic 
study underway, w e sele cted eight clusters from the sample 
of Vikhlinin et al. (1998) with the lowest X-ray fluxes, Lx < 
1.5 X lO'^'^/ig,)^ ergs s^', in the redshift range z — 0.2-0.3. A 
detailed analysis of these eight clusters, including results from 
extensive follow-up spectroscopy, which is in preparation, will 
provide a cruci al link between studies of mas sive clusters (e.g. 
Balogh et al. 1999; Po ggian ti et al. 1999 ) and groups (e.g. 
Zabludoff & Mulchaey [1998| ; Tran et al. [200 1[ ). Comple- 
menting these are nine clusters with high X-ray luminosities 
Lx > ^O'^^hJo ergs s"' at 0.1-2. 4 keV , and taken from the 
XBACS sample of Ebeling et al. ( [1996[ ). The mean Lx of the 
High-Lx clusters is a factor of 30 larger than that of the low 
luminosity cluster sample. It is expected that this difference 
in Lx directly corresponds to a difference in mass, since the 
correlation (Lx oc M^) is well established observationally and 



understood theoretically in th i s L, range (^ewis et al. 1999 
ICooray 1999[; [Wu et al. 1999|; [Xue and Wu 200C|; B abul et al 
2001). While the two samples have very different X-ray prop- 
erties, they cover an overlapping range in local galaxy density. 
We can therefore test directly how galaxy properties depend 
separately on cluster mass and local density. 

To enable the cleanest comparison between the two samples 
we have ensured that the observational data collected are as 
homogeneous as possible in terms of depth, filter and detec- 
tor characteristics. Moreover, since the clusters are all at simi- 
lar redshifts, there is minimal variation in k-corrections, spatial 
resolution, background contamination, field of view, and abso- 
lute magnitude limit between them. 

The data samples and measurements are described in §2. 
In §3 we show the dependence of the morphological mix of 
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Table 1 



Log of the //5r observations of the cluster samples 



Name 


R.A. 


Dec. 


z 


Exposure 


Lx 


(0. 1-2.4 keV) 








(J2000) 




time (ks) 


lo-* 


'^h~~ ergs s^' 






Low— Lx Sample 


















CI 08 18+56 


08 18 58 


+56 54 34 


0.248 


7.2 




0.56 


23.00 


86 


CI 08 19+70 


08 19 23 


+70 54 48 


0.230 


6.9 




0.48 


22.80 


34 


CI 0841+70 


08 41 43 


+70 46 53 


0.240 


6.9 




0.46 


22.90 


39 


CI 0849+37 


08 49 11 


+37 31 25 


0.235 


7.8 




0.73 


22.85 


68 


CI 1309+32 


13 09 56 


+32 22 31 


0.294 


7.8 




0.74 


23.40 


99 


CI 1444+63 


14 44 08 


+63 44 58 


0.297 


7.5 




L46 


23.45 


96 


CI 1701+64 


17 01 46 


+64 21 15 


0.222 


7.5 




0.15 


22.70 


64 


CI 1702+64 


17 02 13 


+64 20 00 


0.242 


7.5 




0.33 


22.90 


61 


High-Lji; Sample 


















A68 


00 36 59 


+09 08 30 


0.255 


7.5 




18.7 


23.00 


102 


A267 


01 52 52 


+01 02 46 


0.230 


7.5 




16.9 


22.80 


90 


A383 


02 48 07 


-03 29 32 


0.185 


7.5 




11.7 


22.30 


73 


A773 


09 17 59 


+5142 23 


0.217 


7.2 




16.0 


22.90 


155 


A 963 


10 17 10 


+39 01 00 


0.206 


7.8 




12.6 


22.50 


96 


A 1763 


13 35 17 


+40 59 58 


0.228 


7.8 




11.8 


22.80 


116 


A 1835 


14 01 02 


+02 51 32 


0.253 


7.5 




48.3 


23.00 


141 


A2218 


16 35 54 


+66 13 00 


0.171 


6.5 




10.9 


22.10 


97 


A2219 


16 40 21 


+46 41 16 


0.228 


14.4 




25.1 


22.77 


110 



the galaxies on cluster X-ray luminosity and projected local 
galaxy density. We discuss the observed correlations in the 
context of models in which the disks or bulges of galaxies 
alone are altered, in §4. Our findings are summarized in §5. 
For all cosmology-dependent calculations we assume A — 0.7, 
J7„, = 0.3 (ACDM) and parametrise the Hubble constant as 
//o =50/z5okms-i Mpc"'. 

2. OBSERVATIONS 
2.1. Data 
2.1.1. Low-Lx Sample 

We selected eight X-ray faint clusters in the northern hemi- 
sphere clusters from the sample identified by Vikhlinin et al. 
( 1998 ) in serendipitous, pointed ROSAT PSPC observations. 
The sample was restricted to a relatively narrow redshift range, 
z = 0.22-0.29 {az/z ^ 0.1) and a mean redshift of z = 0.25, 
to reduce the effects of differential distance modulus and k- 
correction effects on the comparison between the systems. 
The X-ray luminosities of these systems range from 0.15 to 
1.5 X IQ^h-^^ ergs s-i[0. 1-2.4 keV] (Table 1). We compute Lx 
in the 0. 1-2.4 ke V band from the observed fluxes in the 0.5-2.0 
keV band, corrected for galactic Hi absorption and assuming 
a k-correction appropriate for an intra-cluster gas temperature 
equal to that ex pected from the local L, — kT relation ( Allen 
and Fabian 1998; |Mai-kevitch 1998D , using the software pack- 
age XSPEC. We adopt cluster redshifts obtained in the course 
of our own spectroscopic follow-up observations, the results 
of which will be published separately; our red shifts agree well 
with the redshifts published in Vikhlinin et al. ( 1998 ). We refer 
to these clusters as the Low-L^ sample. 

The properties of the Low-Lx clusters are summarized in Ta- 
ble 1, where we list the cluster name (column 1), J2000 coordi- 
nates (2,3), mean redshift from our spectroscopic data (4), ex- 
posure time in ks (5) and Lx [0. 1-2.4 keV] in units of h^^ ergs 
s^' (6). We also list in column 7 the magnitude limit adopted in 
our analysis (see §2.3), and in column 8 the number of galaxies 
brighter than this limit. 

For each cluster, three single orbit exposures with WFPC2 



in the F702W filter were obtained with HST during Cycle 8, 
with exposure times ranging from 2100s to 2600 s per orbit. 
The pointing positions for the three exposures were offset by 
10 Wide Field Camera (WFC) pixels from one another; dur- 
ing the reduction procedure they were aligned and coadded to 
remove cosmic rays and hot pixels. The images were not driz- 
zled or regridded, as this does not accurately preserve the noise 
characteristics of the pixels, and the surface brightness fitting 
software that we use is sensitive to the regridding pattern. Af- 
ter coadding the frames, the image from each WFC chip was 
trimmed to the area of full sensitivity. The three WFC chips 
were not mosaicked, and we have not considered the Planetary 
Camera images in this analysis. The photometry is calibrated 
on the Vega system, with updated zero points taken from the 
current instrument manual. The final images reach a 3-cr point 
source sensitivity of Rn)2 ^ 25.5, and cover a field of 2.5' x 2.5' 
(or Q.%h^Q Mpc at z = 0.25) with an angular resolution of 0.17" 
(^l/i5o'kpc). 



2.1.2. High-Lx Sample 

The comparison sample of very X-ray luminous clusters 
comprises nine clusters with > lO'^^h^^^ erg s^' (0. 1-2.4 
keV) taken from the XBACS sample of Ebeling et al. |1996 . 
With redshifts 0.17 < z < 0.26 (mean of z = 0.21) they lie in 
a narrow redshift range very comparable to that of the Low-Lx 
sample. Redshifts and X-ray luminosities are taken from Ebel- 
ing et al. 
Smith et al 

Struble & Rood 1999). The cluster specifics are given in Table 


Our High-Lx clusters are part of a larger sample of 12 for 
which WFPC2 images were obtained (including some archival 
data) in the course of an HST imaging program focused pri- 
marily on gravitational lensing by intermediate redshift clusters 
(see Smith et al. 2001a, 2001b for first results). The parame- 
ters of the HST imaging are the same as detailed before for the 
Low-Lx sample: each cluster is observed with one HST WFPC 
pointing in the F702W filter for three orbits, with the pointing 



1998 (except the luminosity of A383, taken from 



2001a, and the redshift for Abell 1763, taken from 
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for each orbit shifted by 10 WFC pixels to allow the removal of 
cosmic rays and hot pixels (with the exception of the archival 
data for A2219, which is observed for 6 orbits, and A2218, 
which is shifted only 3 pixels between exposures). The reduc- 
tion procedure is identical to that used for the Low-Lx sample. 



2.1.3. Field Sample 

To provide a reference field sample to compare to the clus- 
ters, and equally importantly to allow us to correct the clusters 
for fore- and background contamination, we have also analysed 
images for eight deep WFPC2 fields se lected from the Medium 
Deep Survey^ (MDS, Ostrander et al. |1998| ; Ratnatunga et al. 
1999). These images were obtained from the MDS archive at 



STScI and are akeady reduced, but coadded without rejection 
of hot pixels. This degrades the cosmetic quality and the num- 
ber of useful pixels, but does not otherwise affect our analysis 
because bad pixels are flagged and rejected from the surface 
brightness fits. 



Table 2 
Log of MDS FIELDS. 



Field 


Exposure 
time (ks) 


Ni,.,,{hu < 22.6) 


u2h91 


28.8 


34 


ucilO 


10.8 


27 


umdOV 


9.6 


20 


umdOa 


8.7 


34 


umdOe 


8.4 


23 


umdOh 


8.2 


28 


ustOO 


12.1 


41 


uwpOO 


8.4 


24 



The MDS images were taken with the F606W and F814W 
(/814) filters. We choose to analyse the F814W images, rather 
than the F606W images, as the k-correction to F702W (used 
for the cluster observations) from F8 14W is typically more uni- 
form over the redshift range of interest [z < 0.8), compared with 
that from F606W, for which the 4000A break moves into the 
F606W passband at z ^ 0.2. To estimate the field contamina- 
tion in our cluster samples we must calculate the surface density 
of field galaxies to a fixed /?702 mag nitude limit. From the deep 
number counts of Metcalfe et al. (2001 ), the median color of 
galaxies brighter than Rjo2 = 23 is {Rj()2 — hu) = 0.4, adopt- 
ing the tra nsformation /xi4 ~ I and Rjm = R — 0.2 (F ukugita 
et al. 1995). This is consiste nt with the median gala xy color 
(y - 7 ) pa 1.1 in the CFR S ( [Crampton et al. 1995| ) and the 
MDS ( [Roche et al. 1996| ), at = 22.5. The color distri- 
bution of galaxies at this depth is broadly distributed, with an 
approximate standard deviation of 0.2 mag; none of the results 
in this paper are sensitive to variations in the (/?702 — /814) color 
within this range. We find a field galaxy density of 19440 deg^^ 
brighter than /^h = 22.5 in the MDS imag es, fu lly consistent 
with the number counts of Metcalfe et al. (2001). We neglect 
any correction in the field density due to lensing by the cluster, 
which we expect to be below the level of field-to-field variance 
(ISmith et al. 2001b|). 



be a large effect (e.g. Saglia et al. 200C), and we do not cor- 
rect for it. Also, selection in redder filters tends to preferen- 
tially select bulge-dominated galaxies; however, the results of 
a m orphological analysis o f the MDS in the V6O6 and /gu fil- 
ters (Ostrander et al. 1998) show this differential effect to be 



small. Finally, the k-corrections for the bulge and disk compo- 
nents will be different across our cluster sample. This differen- 
tial effect is small enough that it can be safely ignored over the 
small redshift range of our sample, but care must be taken when 
quantitatively comparing these results with others at different 
redshifts. A difference of 0.1 magnitudes in the k-correction 
for the bulge, relative to the galaxy as a whole, will result in a 
10 per cent error in the fractional bulge luminosity. 

In Table 2 we list the eight MDS fields used, with their ex- 
posure times in column 2, and the number of galaxies above 
the fiducial limiting magnitude for morphological classification 
(/814 = 22.6, see §2.4) in column 3. Across the eight fields we 
find a total of 231 galaxies brighter than /814 ~ 22.6 (Rjui ~ 
23.0), with a field-to-field standard deviation of 7 galaxies. 

2.2. Source Detection 

Sources were dete cted in the HST images usi ng the SExtrac- 
tor software v.2.1.6 (Bertin and Arnouts 1996). For a source 



to be accepted the signal in at least 10 of its WFC pixels (0.1 
arcsec^) had to be a minimum of 1.5-cr above the background. 
The faintest sources which are reliably detected using these 
criteria have Rjo2 ^ 25-25.5. Our results are insensitive to 
these criteria, as structural parameters can only be reliably de- 
termined for galaxies well above the detection limit. Our re- 
sults do depend, however, on the deblending parameters, as the 
surface brightness fits (see §2.3) use the segmentation image 
produced by SExtractor to identify which pixels belong to the 
galaxy. In the first pass, we used 32 deblending sub-thresholds, 
with a minimum contrast parameter of 7.5 x 10^"*, which de- 
blends well the majority of galaxies. However, in some cases 
this process is too effecient, and identifies, for example, bright 
knots in spiral disks as separate sources. An even more trouble- 
some problem arises in crowded regions, when a smaller galaxy 
is deblended from the flux profile of a brighter galaxy. Often 
SExtractor correctly identifies the separate centroids, but incor- 
rectly associates a large number of pixels from the bright galaxy 
with the fainter In total, about 20 per cent of the sources were 
improperly deblended on the first trial. We therefore repeat 
our measurements on these sources, after interactively choos- 
ing more appropriate deblending parameters. 

2.3. Morphology Measurements 

The visual classification of galaxies onto the revised Hub- 
ble system (or variations on it) has remained a central theme 
of recent morphological studies of distant galaxies with HST 
([Abraham et al. 1996b| ; |Smail et al. l"997| [van den Bergh 1997 



[Couch et al. 1998 ). However, the high resolution and uniform 



There are several possible biases resulting from adopting 
F814W observations for our field correction. The first is that 
galaxies have smoother profiles and more dominant bulge com- 
ponents in the redder filter. However, this is not expected to 

'The Medium Deep Survey catalog is based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. The Medium Deep Survey is funded by 
STScI grant G02684. 



quality data of these HST data has also encouraged the devel- 
opment of automated, machine-based techniques, which allow 
quantitative measures of morph ology related to, b ut distinct 
from, the re vised Hubble system (ISchade et al. 1995|; A braham 
etal. 1996a; pdewahn et al. 1996|; [Naim et al. 1997t M arleau 
and Simard 1998; [Brinchmann et al. 1998f ). For the purposes 
of this paper we have chosen to follow the latter approach and 
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Fig. 1 . — Eighteen galaxies (9 bright and 9 faint) randomly selected from the cluster and field samples. Apparent i?702 magnitude and B/T values are shown in each panel. Each image 
is 6" on a side, and contours are arbitrarily spaced. The top row shows galaxies with B/T< 0.1, which are primarily spirals. Galaxies in the middle row have 0.4 < B/T < 0.6, while the 
bottom row shows galaxies with B/T> 0.8. 



hence analyse our sample using machine-based classifications 
based on fits to the two-dimensional surface brightness profiles 
of galaxies in the WFPC2 images. 

To fit the data we use the IRAF package GIM2d v2.2.0, writ- 
ten by Luc Si mard*^. The details of this program are given in 
Simard et al. (2001). GIM2d fits the sky-subtracted surface 
brightness distribution of each galaxy with up to twelve param- 
eters describing a "bulge" and "disk" compo nent. The algo- 



rithm is similar to that used by Schade et al. (1996c; 1996b), 



except that images are not symmetrized before fitting. Clearly, 
this two component model is a simplified approximation to the 
structure of real galaxies; furthermore, the distinction between 
a "disk" and "bulge" component is made purely on the surface 
brightness profile, and does not necessarily corresp ond to kine- 
matically distinct components (see Simard et al. 2001[ ). For 
simplicity we nevertheless use the terms disk and bulge for 
these components in the following discussion. 

The disk component is approximated by an exponential 
model, with parameters corresponding to the scale length, po- 
sition angle and inclination. For the bulge, we fit a simple de 
Vaucouleurs profile instead of the more general Sersic profile 
because the former generally provides a good fit to bright el- 
liptical gal axies and the bulges of early-type spiral galaxies (d e 
Jong 1994; [Courteau et al 1996|; [Andredakis 1998|). Sersic 



pro files are more a ppropriate for the bulges of late-type galax- 
ies (ie Jong 1994); however, these small bulges will be of low 



signal-to-noise in our images, making an accurate determina- 
tion of the Sersic index difficult. 

GIM2d searc hes th e large parameter space of models using 
the Metropolis ( 1953 ) algorithm, which is inefficient but does 
not easily get trapped in local minima. An important step in this 
algorithm is an initial, broad sampling of the parameter space; 
we populate this space with 300 models. 

We adopt a nominal magnitude limit of Rjo2 = 23.0 at z = 



0.25; galaxies brighter than this value are large enough and 
bright enough to be reliably classified. For ACDM, including 
small (< 0.1 mag) k-corrections, this corresponds to M702 = 
-18.2 + 51og/i5(), appr oximately Mr ^ -18 .0 + 51og/i5o, 4.6 
magnitudes below Mr (Blanton et al. 2001). Accounting for 
the range of redshifts spanned by our clusters, we adjust this 
limit so that the data are matched to the same limiting abso- 
lute magnitude. Therefore, for our most nearby cluster A 22 18 
(z — 0.17) we use a magnitude limit of 7^702 = 22.1, while for 
the most distant cluster CI 1444H-63 (z = 0.29) we use Rjoi = 
23.45. We use a field galaxy sample cut at the corresponding 
/814 limit based on a typical color of {Rj[)2 — /814) ^ 0.4, as dis- 
cussed in §2.1. 

We start from an initial list of sources from the SExtractor 
catalogs which are more than 3" from a WFC chip edge, are 
entirely contained within the trimmed image, and are brighter 
than the relevant magnitude limit for the cluster listed in Ta- 
ble 1 . We then exclude sources with half-light radii less than 
0."15, which are stars, and giant arcs in the High-Lx sample. 
The exclusion of the arcs, which are clearly background galax- 
ies, has no significant impact on our results because the back- 
ground contamination in the High-Lx clusters is very small. 

For each galaxy, we run GIM2d using the pixel membership 
defined by SExtractor's segmentation image. This provides best 
fitting estimates and errors for the bulge and disk components 
of the galaxy. It is quite remarkable how good a fit this simple, 
two component model provides to the majority of the galaxies; 
we find reduced values between one and two for more than 
90 per cent of the galaxies. No restrictions regarding good- 
ness of fit or symmetry are imposed in the subsequent analy- 
sis. From the fitted parameters we calculate B/T, the fraction 
of the luminosity in the bulge, relative to the total luminosity. 
This parameter is known to correlate with Hubble type (Simien 
and de Vaucouleurs 1986), though there is considerable scatter 



http://nenuphar.as.arizona.edu/siniard/gim2d/gim2d.htmI 
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about the mean relation. The uncertainties in B/T are deter- 
mined from Monte Carlo simulation in the GIM2D software, 
and have a median value of 0.07. 

To ensure consistency with the results of the photometric fit- 
ting, the total magnitude for each galaxy is taken to be the total 
flux in the model, as computed by the GIM2D software. These 
total fluxes are systematically brighter than the Kron-type mag- 
nitudes calculated with SExtractor, by ~ 0.15 mag, with a dis- 
persion of ~ 0.25 mag. The offset is larger, ~ 0.25 mag, for the 
brightest galaxies, R-iq2 < 18; fainter than this magnitude there 
is no significant trend with magnitude. 

In Figure 1 we show 6" images of 18 galaxies, randomly 
selected from the Low-L^, High-Lx and MDS samples. The 
RiQ2 magnitude and BAT value are printed in each panel. In the 
top row, three bright and three faint galaxies with B/T< 0.1 are 
shown. This regime of B/T consistently selects disk-dominated 
galaxies, despite the common occurrence of spiral arms and 
bright HII regions. This is not too surprising, since the disk 
decomposition is based on the assumption that the disk surface 
brightness follows an exponential profile, as is the case in the 
majority of spiral galaxies (e.g. de Jong 1996). In the next 
two rows, we show images of galaxies with 0.4 < B/T < 0.6, 
and B/T > 0.8, respectively. This allows the reader to asso- 
ciate B/T measurements with the more traditional morphology 
of the galaxy (see also Marleau & Simard et al. 1998; Tran et 
al. 2001; Simard et al. 2001). 



B/T distributions for the individual Low-Lx clusters are sta- 
tistically consistent with the mean distribution, so there is no 
significant variation from cluster to cluster. 
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Fig. 2. — The luminosity functions of galaxies in the Low-Lx (heavy, solid line) and 
High-Lx (heavy, dotted line) clusters after field coiTection, both renormalized to the av- 
erage number of galaxies per cluster in the High-Lx sample. The error bars include the 
unceitainty in the background correction; the field component subtracted from the Low-Lx 
and High-L;^ samples are shown as the thin solid and thin dotted lines, respectively. 



3. RESULTS AND ANALYSIS 

3.1. Morphological Distributions 

The luminosity functions of the Low-Lx and High-Lx sam- 
ples are shown in Figure 2, normalized to the average number of 
galaxies per cluster above the magnitude limit in the High-Lx 
sample. The background subtracted is also shown in this Fig- 
ure. We find no gross dependence of the luminosity function 
on environment, in agreemen t with the results at lower redshift 
from Zabludoff & Mulchaey (2000). We can therefore compare 
the distributions of B/T in the High-Lx and Low-Lx clusters, 
including all galaxies brighter than M702 = — 18.2 + 51og/i5o 
(> O.OIL*). Since the redshift distributions of the clusters (Ta- 
ble 1) and the luminosity functions of their constituent galaxies 
are very similar, any differences seen in these distributions must 
therefor e be related (though not necessarily directly, as we dis- 
cuss in §3.2) to the difference in X-ray luminosity, cluster mass 



or local galaxy density. 

To correct for the fore- and background contamination of the 
clusters, we first consider the distribution of B/T in the field 
galaxy sample from the eight MDS fields, shown in the top 
left panel of Figure 3. The distribution has a strong peak at 
B/T= 0, and 65 ±9 per cent (error from the field-to-field stan- 
dard deviation) of all galaxies have B/T< 0.4; i.e., they are 
disk-dominated. In the remaining panels on the left side of Fig- 
ure 3 we show the distribution of B/T in each of the eight Low- 
Lx clusters (dotted histograms). We correct for foreground and 
background contamination by subtracting the mean field B/T 
distribution; the corrected distribution is shown as the solid his- 
tograms. In most cases, the effect of background contamination 
is small, and generally serves to reduce the fraction of disk- 
dominated galaxies by only a few per cent. In most clusters the 
field-corrected B/T distribution is reminiscent of that seen for 
the field population, with a comparatively flat distribution ex- 
cept for a pronounced peak in several clusters at B/T ^ 0. The 



The B/T distributions for the clusters of the High-Lx sample 
are shown in the right hand panels of Figure 3, again with and 
without field correction (this is less important for this sample 
than for the Low-Lx clusters). In general, the High-Lx sam- 
ple shows more marked differences from the field population 
than the Low-Lx systems. In particular, most of the clusters 
lack large numbers of galaxies with B/T ~ and instead show 
a broad peak of bulge-strong galaxies with ^IT ^ 0.5. There 
are three notable exceptions to this general statement: A 963, 
A 267 and A 383, which all show a peak at B/T=0, and distri- 
butions which are broadly similar to the Low-Lx clusters. We 
discuss these three discrepant clusters in more detail in §Q 

The simplest way to characterise the morphological mix of 
the clusters is to evaluate the fraction of disk-dominated galax- 
ies in each cluster after field-correction. We arbitrarily define a 
disk-dominated galaxy to have B/T< 0.4; however, our qualita- 
tive conclusions are unchanged for any definition ranging from 
B/T< 0.1 to B/T< 0.5. In Figure 4 we show the fractions of 
disk-dominated galaxies in the clusters as a function of their 
X-ray luminosities. The disk-dominated galaxy fraction in the 
field sample has a mean of 0.65 and we show this for reference; 
along with the 1-cr field-to-field standard deviation of 0.09. The 
error-weighted mean fraction of disk-dominated galaxies in the 
High-Lx clusters is 0.23±0.10, where the quoted uncertainty 
represents the \-tj standard deviation from cluster to cluster. 
This 10 per cent variation is well above the statistical uncer- 
tainty, and is due to the three clusters noted earlier (A 963, 
A 267 and A 383) which have unusually high fractions of disk- 
dominated galaxies. For the Low-Lx cluster sample the spread 
in B/T is similar, but the disk-dominated fraction is much higher 
on average, at 0.40±0.09. However, even in these systems the 
disk fraction is still significantly lower than the general field. 
This suggests that the processes influencing galaxy morphol- 
ogy can effect galaxies inhabiting structures less massive than 
X-ray luminous clusters, in agreement with the results of Tran 
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Fig. 3. — Left: The distribution of B/T in the eight Low-Ly clusters, and the mean distribution of the MDS field sample. For the eight clusters, the dotted histograms show the 
distribution of all galaxies, and the solid histograms are corrected for fore- and background contamination by subtracting the field distribution. Right: The B/T distributions of the High-Lx 
clusters. 



et al. (2001). However, we caution that the field galaxy sample 



is magnitude limited, and therefore has absolute magnitude and 
redshift distributions different from those of the volume-limited 
cluster samples. In particular, we expect the redshift distribu- 
tion of the field sample to peak ai z~ 0.55 ( Crampton et al. 
1995). 



d 

V 
E- 

o 



ID 
O 



O 

CV! 

Cd 
O 

o 

(0 

u 
6- 



—\ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

Field 



_L 



'■h 



_L 



J I I I I L 



43 44 45 

log [0.5-2.5 keV] (hg^ ergs/s) 

Fig. 4. — The fraction of disk-dominated galaxies (those with B/T<0.4) in the Low- 
Lx and High-Lx samples as a function of the cluster X-ray luminosities {open circles). 
Error bars on the data points are \-a bootstrap estimates, and include the uncertainty in 
the background coirection. The solid diamonds show the eiTor-weighted mean fractions 
for the combined Low-Lx and High-Lx samples, with an eiTor bar which represents the 
1 — cr standard deviation in this fraction from cluster to cluster {i.e. it is not the error on the 
mean). The solid line is the coiTesponding fraction in the field population, and the dashed 
lines show the l-cr field-to-field standard deviation of this quantity. 



The results of the previous section show that the distribution 
of galaxy morphologies is different in the Low-L^ and High- 
Lx cluster samples. We now wish to test whether this is due 
to differences in the local galaxy environment, or whether it 
instead reflects the influence of the global environment on the 
galaxy populations. 

The clusters in the two samples we analyse have significant 
differences in their physical properties, with the mean Lx of 
the High-Lx sample a factor of ^ 30 times higher than that 
of the Low-Lx sample. This corresponds to a difference of a 
factor 5 in mass, or a factor of 2 in vir ial radius (from 
/J^jj.; e.g. see Babul et al. [200 1| ). As the galaxy 
samples in the High-Lx and Low-Lx clusters are selected from 
the same physical region (~ 3' in diameter or ^ I/i^q' Mpc for 
a ACDM cosmology), the galaxies in the Low-Lx sample are 
selected from a region covering a substantially larger fraction 
of the virial radius of the clusters, compared with the High-Lx 
sample. Within relaxed systems local density is strongly corre- 
lated wth^admsjiorm^li^ed^ et al. 
1997); we therefore expect that the typical galaxy density in 
the regions of the High-Lx and Low-Lx clusters analysed here 
will differ. Thu s, if a strong m o rphology-density relation exists 
in the clus ters ([Dressier 1980|; Postman and Geller 1984|; T ran 
et al. 2001), the systematic variation in the local galaxy density 
between the two samples will lead to differences in the apparent 
morphological mix. 



3.2. Local or global environment? 
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Fig. 5. — Top panel: The distribution of projected local galaxy surface density in the 
Low-Lx (solid line) and High-Lx (doited line) samples. Bottom panel: The fraction of 
galaxies with B/T<0.4 as a function of projected local galaxy density. The values are 
corrected for background contamination as described in the text. The field fraction (which 
corresponds to Siocal ~ 0) is shown as the star at an arbitraiy density for display purposes. 
The circles and triangles correspond to the Low-Lx and High-L^ samples, respectively. 
Bin sizes are adaptively varied so that each includes 75 galaxies; the horizontal error bar 
shows the bin span, and the point represents the median. Vertical eiTor bars are deter- 
mined from bootstrap resampling which accounts for the uncertainty in the background 
correction. 



The local projected galaxy density around eac h gal axy is 
computed in a manner similar to that of Dressier (198C). The 
five nearest neighbor galaxies are found, down to the appropri- 
ate limiting magnitude (see §2.4) and the encompassing area 
is taken to be a circle which extends out to the most distant 
neighbor (correctly accounting for the smaller area available for 
galaxies near the boundary of the WFC mosaic). We choose the 
fifth nearest neighbor, rather than the more traditional choice 
of the tenth, since the latter choice results in an uncomfort- 
ably large smoothing scale, relative to the small area of the 
WFC images. The average field density computed at the ap- 
propriate magnitude limit (S/ = 5.4 arcmin^^ at /814 — 22.5) 
is subtracted to give the local density Eiocai, which is converted 
to physical units {h^Q Mpc^^) assuming a ACDM cosmology. 
When computing the projected densities, all of the surrounding 
galaxies in the photometric catalog above the magnitude limit 
are considered, though we have only determined morphologies 
for galaxies which are well clear of chip boundaries. We have 
not accounted for the fact that the trimmed chip images do not 
exactly join, so local densities near the common boundary of 
two WFC chips will be slightly underestimated. 

The contribution of the background to the observed B/T dis- 
tribution depends on the local projected density of galaxies un- 
der consideration; in denser regions, the background correc- 
tion will be proportionally smaller For the total cluster sam- 
ple, we computed the mean projected field density from the 
MDS, Sy, to the appropriate magnitude limit for each clus- 
ter, and subtracted the corresponding number of galaxies. For 
galaxies within a restricted range of projected densities, we 
scale the MDS B/T distribution appropriately so that a fraction 
/ = S// (Sy + S) of the galaxy population is statistically sub- 
tracted from the observed B/T distribution, where E is the me- 



dian,_ background-corrected density of the galaxy subsample. 
For E ~ 25/1^0 Mpc^^ / ^ 60 per cent, while for E - 250/j^o 
Mpc~^, / 10 per cent. 

In the top panel of Figure 5 we show the Eiocai distributions 
of the Low-Lx and High-Lx samples. As expected, most of 
the galaxies in the Low-Lx sample are drawn from regions of 
much lower galaxy density than seen for the bulk of the galax- 
ies in the High-Lx sample. Note that the typical densities in the 
High-Lx sample are much larger than those shown in Dressier 
( |1980| ), since we are including galaxies about 2.5 magnitudes 



less luminous than in Dressler's study. 
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Fig. 6. — The distribution of B/T ratios is shown for low (lowest 5-25 percentile) and 
high (highest 5-25 percentile) density regions separately, as indicated. The Low-Lx clus- 
ter data are shown in the top panel, and the High-Lx cluster data are shown in the bottom 
panel. Uncertainties are 1 — cr and include the error in background subtraction. 



The morphology-density relation is shown in the bottom 
panel of Figure 5. We use an adaptive bin size so that each bin 
contains 75 galaxies, and compute the field-corrected fraction 
of those that are disk-dominated, B/T< 0.4 (as this definition of 
disk-dominated is arbitrary we will consider the full B/T distri- 
bution as a function of density below). In the individual cluster 
samples there is only weak evidence for a morphology-density 
relation due to the small samples and limited dynamic range 
in density. In Figure 6 we compare the background-subtracted 
distribution of B/T in low (5-25 percentile) and high (75-95 
percentile) density regions within the individual Low-Lx and 
High-Lx samples. The corresponding density limits in each 
sample are also shown in Figure 6. This compares galaxies in 
well separated density ranges above and below the mean of the 
sample, while avoiding the extreme tails of the distribution. No 
statistically significant difference can be seen between the low 
and high density regimes, in either cluster sample. There is 
an indication, in the High-Lx clusters, that the lowest density 
galaxies are distributed toward lower B/T, but the reduced 
statistic for the difference between these two distributions, ac- 
counting for the uncertainty in background subtraction, is only 
1.2''. The T-E is clearly weak, at best; this can be attributed 
partly to the limited dynamic range of densities explored; HST 
imaging over a wider field would allow the extension of the 



analysis of the High-Lx clusters to much lower densities, and 

'We are forced to bin the data in order to perform the background subtraction, and therefore we use the x^ statistic, rather than the Kolmogorov-Smimov statistic. 
The disadvantage is that this statistic is quite sensitive to the chosen binning, and we do not, therefore, consider this a very robust estimate of the significance of our 
result. This is true of all estimated likelihoods based on the statistic, throughout this paper. 
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improve the comparison. 

There is some evidence from Figure 5 that, at a given 
overdensity, disk-dominated galaxies are more common in the 
Low-Lx clusters. We investigate this possibility by restricting 
the galaxy samples in both High-Lx and Low-Lx clusters to 
a common density range where both samples have good statis- 
tics, Eiocai = 50-200. In this density range, there is good over- 
lap between the two datasets, and limited gradient in morpho- 
logical composition. Although this density range is still fairly 
large, our results do not qualitatively change as this interval is 
reduced, at the necessary expense of lower statistical signifi- 
cance. 

The B/T distributions in the High-Lx and Low-Lx clusters, 
restricted to Eiocai = 50-200, are shown in Figure 7. The two 
distributions appear to differ, with a relative excess of galax- 
ies in the High-Lx sample with B/T~0.5, and a corresponding 
deficit of those with BAT ^ 0. A statistical comparison of the 
two normalised distributions gives a reduced of 3.1 and a 
corresponding likelihood of 98%. This difference is much more 
pronounced and significant (99.99% likelihood) if w e exclude 
the three unusual High-Lx clusters discussed in §3.1 ; however, 
we are reluctant to do so without identifying the origin of the 
difference between these and the majority of the High-Lx sys- 
tems. We also reemphasize that the x- test is sensitive to the 
binning of the data and thus does not provide a complete statis- 
tical description of the difference between the samples. 



4. DISCUSSION 

From Figures 3 and 4, it is evident that all of the Low-Lx 
clusters have similar galaxy morphology distributions; all eight 
clusters are dominated by galaxies with BATw 0, and are statis- 
tically indistinguishible from the mean B/T distribution of the 
combined sample. Most of the High-Lx clusters differ from 
this distribution, as they are dominated by galaxies with inter- 
mediate B/T. The exceptions are three clusters, A 963, A 267 
and A 383, which have B/T distributions more similar to those 
found in the Low-Lx clusters. There is nothing strikingly un- 
usual about the morphology of these three clusters; all are dom- 
inated by a single cD galaxy and do not show obvious signs of 
irregularity. The redshifts of these three clusters are not at either 
extreme of the distribution, nor are the galaxy luminosity func- 
tions, or distributions of local projected density, significantly 
different from that of the average High-Lx cluster. 

One possible explanation for the difference could be that 
these three clusters are less massive than the other High-Lx 
clusters, and that their X-ray luminosities are unusually high, or 
overestimated. Alternatively, if the clusters have an unusually 
large number density of fore- and background galaxies along 
the line of sight, the background correction will have been un- 
derestimated. However, the local galaxy density distribution of 
the three clusters is indistinguishible from the High-Lx clus- 
ter average; that is, the galaxy density is as high as we expect 
for the cluster's luminosity. It would seem, therefore, that both 
the luminosity and the density would have to be artificially en- 
hanced for these explanations to work; this would be an uncom- 
fortable coincidence. It therefore appears that galaxy morphol- 
ogy is sensitive to factors other than just local density or X-ray 
luminosity, though it is not clear what they are. 




Fig. 7.— The B/T distribution for galaxies with Sioci = 50-200 Mpc"^. The solid 
and dotted histograms correspond to the Low-Ljt and High-Lx samples, respectively. 



Putting aside for the moment the intrinsic variation in the 
High-Lx sample, we can explore differences in the average 
properties of the High-Lx and Low-Lx clusters. There is some 
evidence that disk-dominated galaxies are relatively more com- 
mon in the Low-Lx clusters, at fixed overdensity, relative to the 
High-Lx clusters (Figure 7). Although this result is of marginal 
statistical significance, it is still of interest to investigate which 
physical process could give rise to the effect. In currently pop- 
ular models of galaxy formation, bulge-dominated galaxies are 
generally created in one of two ways. The first is by the merger 
of two galaxies, in which some or all of the baryonic mate- 
rial (from the progenitor bulges and disks) is converted into a 
new bulge component. In this case, differences in the bulge lu- 
minosity function will reflect differences in the galaxy merger 
histories. Alternatively, a bulge-dominated galaxy can be cre- 
ated by reducing the stellar mass of the disk component, either 
by destroying it (for example, through harassment) or by caus- 
ing it to fade sufficiently following the cessation of star forma- 
tion (as might be expected if ram pressure stripping removes 
all the gas from the disk). In this case, the bulge luminosity 
function should be unaffected by the efficiency of the destruc- 
tion process, but the B/T measure will increase as the disk fades 
or is destroyed. We will attempt to distinguish between these 
two possibilities by examining the relationship between bulge 
luminosity and B/T in the two cluster samples. 

We show the bulge and disk luminosity functions separately 
in Figure 8, for those galaxies from Low-Lx and High-Lx sam- 
ples lying within the density range Eiocai — 50-200. The lu- 
minosity functions are normalized to the average number of 
galaxies in this density range per cluster in the High-Lx sam- 
ple. While the disk luminosity functions appear to be very sim- 
ilar in both cluster samples, there is an indication (not statis- 
tically significant, however) that the bulge luminosity function 
is steeper in the High-Lx clusters; i.e., there is a higher pro- 
portion of low-luminosity bulges compared with the Low-Lx 
sample. This suggests that the difference in B/T distributions 
is reflecting a difference in the properties of the bulges, rather 
than the disks. With some effort, this can be seen more clearly 
in the correlation of B/T with disk or bulge luminosity, shown 
in Figure 9. We again restrict the galaxy population to those 
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Fig. 8. — Left: The luminosity function of tiie bulge component, in the Low-Lx sample (heavy, solid line) and the High-Lx sample (heavy, dotted line), restricted to regions of 
overdensity 50 < Siocai < 200. 0{M) is normalised to the average nuinber of galaxies per cluster in the High-Lx sample, restricted to this density range. The thin solid and dotted lines 
show the subtracted background component. Right: The same, but for the disk components. 



with Siocai ~ 50-200, to minimize morphology-density effects. 
Concentrating first on the top panels, we see the correlation be- 
tween B/T and bulge luminosity, for the Low-Lx and High-Lx 
samples separately. Galaxies with B/T< 0.1 have been arbitrar- 
ily plotted at BAr= 0.1 to make them visible in this logarithmic 
plot. Note also that the shape of the sample luminosity limit, 
combined with the shape of the luminosity function, gives rise 
to some of the structure seen in this figure. It is clear (as it is in 
Figure 7) that the Low-Lx clusters are dominated by galaxies 
with low B/T, while this is much less the case in the High- 
Lx clusters. We can now explore the two scenarios described 
above. If galaxy disks are preferentially stripped, or otherwise 
fade, in the High-Lx clusters, the galaxy density will be shifted 
directly upward in this figure, relative to the Low-Lx distribu- 
tion (see the aiTows labelled "S"). This would serve to push 
many of the galaxies with B/T< 0.2, which are mostly faint, 
below our luminosity limit. In this respect, such a mechanism 
could provide a viable explanation of the difference in galaxy 
distributions seen in these top figures. However, if we now look 
at the bottom panels, it can be seen that this explanation is un- 
likely. In these figures, we plot disk luminosity as the x-axis. 
Here, it is clear that the faint end of the disk luminosity func- 
tion is primarily defined by galaxies with B/T< 0.2. If these 
disks were to fade beyond the luminosity limit, the disk lumi- 
nosity function would necessarily become much shallower, and 
this does not appear to be the case (Figure 8). Instead, the faint 
end of the disk luminosity function in the High-Lx clusters is 
relatively more populated by galaxies with BAT^ 0.5 and disk 
luminosities roughly equivalent to the disk luminosities of the 
B/T^ population in the Low-Lx clusters. 

On the other hand, if the difference between the Low-Lx and 
High-Lx clusters is due to the fact that the bulges are system- 
atically brighter in the High-Lx clusters, due perhaps to a more 
extensive merger history, galaxies in these panels will move di- 
agonally, as indicated by the arrows labelled "M". This expla- 



nation appears consistent with the observations. The B/T< 0.2 
population which dominates the Low-Lx clusters can be trans- 
lated in the sense of brightening their bulges at constant disk 
luminosities, to reproduce the distribution seen in the High-Lx 
clusters. 

We conclude that is unlikely the B/T^^ 0.5 galaxies in the 
High-Lx clusters have been formed from a process which op- 
erates uniquely in this environment by stripping or otherwise 
destroying the disk in a lower B/T system. Rather, the results 
seem to suggest that there are differences in bulge growth in 
the Low-Lx and High-Lx samples. Thi s prov ides quantitative 
support for the suggestion by Dressier ( 1980 ), that the early- 
type galaxies in clusters cannot all be formed by destroying the 
disks of a normal population of late-type galaxies. However, 
this is not evidence that ram pressure stripping (or other pro- 
cesses with similar effect) do not take place at all; at least, they 
could play a strong role in establishing the local morphology- 
density relation. But the difference between the T-S for Low- 
Lx and High-Lx clusters indicates that there are other processes 
at work. 



5. CONCLUSIONS 

We have presented an HST-hased morphological analysis of 
galaxies in clusters from two samples at z ^ 0.25. The two data 
sets are extremely well matched in the properties and quality 
of their data; all observations are 3-orbit HST exposures in the 
F702W filter, and both cluster samples have a similar mean red- 
shift and small redshift range, thus minimizing differences in k- 
corrections and absolute magnitude limits. The only significant 
difference between the samples is the X-ray luminosity of the 
clusters contained in them; the clusters in the Low-Lx sam- 
ple have Lx(0.1-2.4keV)=0. 15-1. 5 X 10^4 ergs s-i, while the 
High-Lx clusters have Lx(0.1-2.4keV)> 10"*^ ergs s^'. These 
data are therefore ideal for testing the dependence of galaxy 
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Fig. 9. — Top: The correlation between B/T and bulge luminosity is shown for galaxies in the High-Lx clusters (left) and the Low-Lx clusters (right), after background subtraction. 

The sample is restricted to those galaxies in the local projected density range Siocal — 50-200 Z?^,, Mpc^^. The density of galaxies in this plane is represented as a contour plot; the contours 
are equally spaced in intervals of 5 galaxies per cell (starting with 1), where a cell is 1 mag and 0. 1 in B/T. Galaxies with B/T< 0.1 are arbitrarily fixed at B/T-0. 1 to be visible in these 
plots. The short-dashed lines show the bulge luminosity of an M* galaxy, and a galaxy three magnitudes fainter, as a function of B/T, for reference. The long-dashed line represents the 
luminosity limit of the sample. The arrows labelled "S" and "M" show the magnitude and direction a galaxy with B/T=0.1 woixld move, if its disk faded by one magnitude (stripping: S) or 
its bulge brightened by one magnitude (merging: M). Bottom: Same as the top panel, but for the correlation between B/T and disk liuninosity. Galaxies with B/T> 0.9 are arbitrarily fixed 
at B/T= 0.9 to be visible in these plots. 
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morphology on local density and the global environment. We 
summarize our findings as follows: 

• Within a fixed physical region (/^ ^hJo Mpc diameter), 
the fraction of disk-dominated galaxies is strongly de- 
pendent on Lx- In the High-Lx sample, the fraction of 
galaxies with B/T< 0.4 is 0.25±0.10, while in the Low- 
Lx sample it is 0.44±0.07. 

• In both cluster samples, the fraction of galaxies with 
B/T< 0.4 is significantly lower than the fraction of such 
galaxies in the field, 0.65±0.09 (where the uncertainty 
represents the field-to-field variation). The origin of this 
discrepancy is, however, unclear because the comparison 
is between a volume-limited sample at z 0.25 (the clus- 
ters) and a magnitude-limited one with a broad redshift 
distribution peaked at z 0.55 (the field). 

• There is a dependence of the typical B/T of a galaxy on 
its local projected galaxy density. This is partly responsi- 
ble for the above result, as galaxies in the High-L^ clus- 
ters are generally located in denser environments. How- 
ever, when we compare galaxies at similar local densities 
in the Low-L^ and High-L^ samples, small differences 
persist, with the High-L^ clusters showing an excess of 
bulge-strong galaxies compared to the Low-L^ sample. 
This difference is moderately significant, with confidence 
limits of 98% (as determined by a test, which is sen- 
sitive to how the data are binned). 

• To investigate this behaviour in more detail we focus on 
the bulge luminosity function in the two samples. We find 
some evidence for a difference in the bulge luminosity 
function of galaxies between the Low-Lx and High-Lx 
clusters when we restrict the comparison to galaxies with 
local densities between 50 and 200 /150 Mpc^^, although 
the difference is not highly statistically significant. 



• Three of the High-Lx clusters show B/T distributions 
which are more similar to the average Low-Lx B/T dis- 
tribution than the majority of High-Lx clusters. Galaxy 
morphology therefore appears to depend on another, un- 
known parameter, besides local density and X-ray lumi- 
nosity. 

These observations suggest that the factors influencing 
galaxy morphlogy are: 1) efficient in relatively low-mass clus- 
ters; and 2) influence the luminosity of the bulge component, 
and not the disk alone. The most likely interpretation in our 
opinion is that the merger histories of galaxies in low and high 
mass clusters are different, in that merging has played a larger 
role for galaxies 1-3 magnitudes below ^ L* in more massive 
clusters. This mig ht be expected from models of hierarchi - 
cal cluster growth ( Lacey and Cole 1994 ; Kauffmann 1996 ), 
since the merger history of massive clusters tends to be more 
extended. On the other hand, galaxy-galaxy mergers are un- 
common in regions of high velocity dispersion, and it is un- 
clear what variation in the bulge luminosity function would be 
predicted by the models. 
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